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Abstract 

In frame of a phenomenological approach based on Compton-like Feynman amplitudes, we study 
the annihilation channel in antiproton nucleon collisions with production of a single charged or 
neutral pion and a lepton-antilepton pair. These reactions allow to access nucleon and axial elec- 
tromagnetic form factors in time-like region and offer a unique possibility to study the kinematical 
region below two nucleon threshold. The differential cross section in an experimental set-up where 
the pion is fully detected is given with explicit dependence on the relevant nucleon form factors. 
The possibility to measure heavy charged pion in the annihilation channel is also discussed. 
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I. INTRODUCTION 



In this paper we study the annihilation reaction p + n — > tt~ + £~ + £ + and p + p — > 
7r° + £~ + £ + , £ — e, fM which is the crossed process of pion electroproduction on a nucleon 
N: e~ + N-^e~ + N + 7T. It contains the same information on the nucleon from factors 
in a different kinematical region. This process is also related to the pion scattering process 
ii + N—>-N + £~+£ + which was first studied in pj. In this work it was already pointed 
out that the p + p annihilation process with pion production, under study here, renders 
possible the determination of the nucleon electromagnetic form factors in the unphysical 
region, which is otherwise unreachable in the annihilation reaction p + p — > e + + e~. In 
Q] a general expression for the cross section was derived and numerical estimations were 
given in the kinematical region near threshold. In this paper we extend the formalism in 
two directions. We take into account a larger set of diagrams which can contribute and give 
special emphasis to the possibility to access the axial nucleon form factors in the time-like 
region. For the annihilation process pN — > £ + + £~ + 7r°'~, with £ = e, or /i, we will consider 
two reactions: 

P(Pl) +P(P2) -> TT ^) +i + {p+) +t~(P-), (1) 

and 

p(pi) + n(pa) -> 7r-(g w ) + £ + (p+) + l~ (p_), (2) 

where the notations of the particle four momenta are indicated in brackets. The interest 
of these processes lies in the possibility to access nucleon and axial form factors (FFs) in 
the time-like region. The expected values of the total cross sections at energies of a few 
GeV drop with the lepton pair invariant mass q 2 , but are of the order of several nb and up 
to mb below threshold. Therefore these reactions are especially interesting as they will be 
measurable in next future at hadron colliders (or the crossed reactions at lepton colliders). 
The present work aims to evaluate the differential cross section for experimental conditions 
achievable at the future FAIR facility 

A. Electromagnetic (Vector) Form Factor 

The structure of proton and nucleon is of great interest since it explores the quantum 
field theory of strong interaction, quantum chromo dynamises (QCD), in a region where 
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the interaction between the constituents of quarks and gluons can not be treated as a 
perturbation. Its properties have been explored in many observables. For example, the 
electromagnetic form factors of proton and nucleon are measured in elastic electron-nucleon 
scattering with a space-like four momentum transfer q 2 < via the process e + p — > e + p 
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. For a momentum transfer of —1 GeV 2 < q 2 < GeV 2 , the electromagnetic form factors 
are well known for both neutron and proton 5|. 

The interest in the form factors of the nucleon for the region of q 2 < — 1 GeV 2 has recently 
been renewed by a measurement of the ratio of electric form factor over the magnetic from 
factor [iJGeI Gm a t TJNAF in a g 2 -range up to - 5.8 GeV 2 . The method of polarization 
transfer [6j was employed for the first time at high values of q 2 using a longitudinally polarized 
electron beam and measuring the polarization of the recoil proton |7|. The experiments 
gave the surprising new result that the ratio \x v GeI Gm decreases from unity substantially 
for high q 2 - values exhibiting eventually even a zero-crossing around q 2 ~ —8 GeV 2 . This 
new result is in contradiction with all data obtained from a Rosenbluth separation fit of 
unpolarized measurements |8] giving a ratio of h p Ge/Gm ~ 1- The present understanding 
of this discrepancy is that the analysis of the unpolarized data lacks important corrections 
stemming from neglected higher order radiative corrections jsj or electromagnetic processes 
like two-photon-exchange [10]. 

The situation of the experimental determination of the electromagnetic form factors in 
the time-like domain (i.e. q 2 > 0) is quite different. Existing data have explored the two 
basic processes: e~ + e + — > p + p and p + p — > e~ + e + . The precision of the cross section 
data in the time-like region is much lower and a determination of the individual electric and 
magnetic form factors has not been done with sufficient precision. 

For both processes there is a threshold in q 2 for producing the two nucleons at rest or 
annihilation of the two nucleons at rest which amounts to q 2 > (m p + nip) 2 = 4m 2 . In 
the region of < q 2 < 4m 2 (sometimes called the "unphysical region") there are no data 
available. On the other hand this g 2 -interval is of great interest since the intermediate virtual 
photon as well as the pp-pair can couple to vector meson- and pp-resonances respectively 
and thereby enhance the form factors substantially This mechanism has been 

used to explain the large cross section in the time-like region above the threshold q 2 > 4m 2 . 
Another possibility to explain this large cross section above threshold has been used recently 
by showing, that the large cross section in the annihilation process p + p — > e + + e~ can be 
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related to the pp scattering length [3] . 

The work presented here is focused on a possibility of extracting the electromagnetic 
form factors of the nucleon in the "unphysical region" and above using the process p + p — > 
it + e + + e - , as it has been proposed earlier 0, [J. The measurement of the cross section of 
this reaction will be accessible at the future FAIR facility at GSI. 

B. Axial Form Factor 

In addition, we explore here also the new idea of accessing the axial vector current of the 
nucleon. The axial form factors in the space-like region are measured in nuclear /3-decay, in 
neutrino scattering, in muon capture and in pion-electroproduction (e + p — > e + n + tt + ). 
A review on the present status of the axial structure of the nucleon in the space-like region 
is given in Ref. [14] . The first methods represent a rather direct measurement in the sense 
that the axial coupling of the weak charged currents is used to measure the axial form 
factor. The extraction of the axial form factor from pion-electroproduction in the space 
like region is possible due to the application of a chiral Ward identity referred to as the 
Adler-Gilman relation 1151. A recent review on the theoretical development can be found 
. Corrections to order 0(p 4 ) have been calculated in the framework of chiral 
perturbation theory in Ref. ^|. There are no data available on the axial form factor in 
the time-like region. A direct measurement would be possible by studying the weak neutral 
or charged current in the annihilation of pp. Such a measurement is not accessible with 
present day experimental techniques. On the other hand an application of the Adler-Gilman 
relation to the matrix element of the crossed channel of pion-electroproduction, namely 
p + n — > 7r~ + e + + e~, renders the possibility of accessing the axial current also in the time- 
like region around and below threshold. Despite the theoretical uncertainties, namely on 
one hand the applicability range of the Adler-Gilman relation which is strict only at the tt- 
threshold and on the other hand the treatment of the off-shell nucleon between the it- vertex 
and the virtual photon vertex, the measurement of the cross section in p + n — > 7r~ + e + + e~ 
at low energies would render the first estimation of the axial form factor in the time-like 
region. Also here, the measurement of the cross section of this reaction will be accessible at 
the future FAIR facility at GSI. 
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C. Approach 

Our approach is based on Compton-type annihilation Feynman amplitudes and aims 
to establish the matrix element of the processes (0) and The main uncertainty in our 
description in terms of Green functions of mesons and nucleons (and their expected partners) 
is related to the model dependent description of hadron FFs and to the modeling of excited 
hadronic states. 

The paper is organized as follows: the formalism is developed in Section II and some of 
the kinematical constraints for the considered reactions are discussed in Section III. Section 
IV contains a discussion on nucleon FFs and of our choices of parametrization, for electro- 
magnetic as well as for axial FFs. In Section V the numerical results for the differential 
cross section of the considered processes will be presented. In Conclusions we discuss the 
results and summarize the perspectives opened by the experimental study of these reactions, 
including possible manifestation of heavy (radial) excited n states. 

II. FORMALISM 

Let us consider the reactions (JTJ) and (J2J and calculate the corresponding matrix elements 
in the framework of a phenomenological approach based on Compton-like Feynman ampli- 
tudes. The Feynman diagrams for the reaction (0), are shown in Fig. and^i, for pair 
emission from the proton and from the antiproton, respectively. For the reaction (J2J), the 
corresponding Feynman diagrams are shown in Figs. EK, Eb an d Et for pair emission from 
the charged pion, from the antiproton and from the neutron, respectively. 

As previously discussed in Ref. Q|, no free nucleons are involved in the electromagnetic 
vertexes of pions and nucleons, as one of the hadrons is virtual, and rigorously speaking, the 
form factors involved should be modified taking into account off mass shell effects. However 
we will use the expression for electromagnetic current involving on mass shell hadrons. A 
discussion of the errors and the consequences of such approximation is given below. 

The vertexes ^y*NN* — > 7*AA and 7*7T7r — ► 7*7T7r* contain the information on the 
electromagnetic form factors (FFs) of proton, neutron and pion (FFs of antiproton differ by 
sign from proton FF, due to charge symmetry requirements). Nucleon FFs can be expressed 
in terms of Dirac and Pauli FFs Ff '^{q 2 ), which enter in the expression of the electromagnetic 
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FIG. 1: Feynman diagrams for the reaction p + p — > tt° + £ + + £ . 



current: 

<iV(p')|r,(g) 7V |iV(p) >=u(p') 



u{p), N — n,p, (3) 



where M is the nucleon mass and q is the four-momentum of the virtual photon. The nucleon 
FFs in the kinematical region of interest for the present work are largely unexplored. The 
assumptions and the parametrizations used for FFs in the numerical applications are detailed 
below. 

The pion electromagnetic FF, F n (q 2 ) is also introduced in the standard way. The corre- 
sponding current has the form: 



(4) 



with qi and q2 ingoing and outgoing charged pion momenta, and q n = q\ — q%. Special 
attention must be devoted to the pion nucleon interaction, in the vertexes ttNN that are 
parametrized as: 



v{pi)l- > u{p 2 )g- K NN{s), and v(p 1 - q)^u{p 2 )g 1T NN{'ml). 



(5) 



with s = (pi + p 2 ) 2 . The vertex of irNN interaction is related to the general axial vector 
current matrix element: 



< N(j/)\A1(0)\N{p) >=u(p' 



G^hp+£GW) + i^G T (q 2 ) 



T; 



J5~^u(p), (6) 



where = p'^ — p M , (^(g 2 ) is the axial nucleon FF, Gp(q 2 ) the induced pseudoscalar FF, 
and Gxiq 2 ) the induced pseudotensor FF. In the chiral limit, the requirement of conservation 
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FIG. 2: Feynman diagrams for the reaction p + n — > ir + £ + + . 



of the axial current leads to the relation: 



AMG A {q 2 )+q 2 G P {q 2 ) = 0, 



(7) 



which shows that Gp(q 2 ) has a pole at small q 2 . Indeed, assuming that the axial current 
interacts with the nucleon through the conversion to pion interaction, one obtains: 

4,Mf^ NN (q 2 ) 



G P (q 2 ) ■ 

q- 

Comparing Eqs. ((Zj) and (JHJ) one obtains the Golberger-Treiman relation: 

G_a _ 9-kNN 

We suggest a generalization of this relation in the form: 

MG A (s 



(8) 



(9) 



9{S = 9nNN{S) 



G A (0) = 1.2673 ±0.0035. 



(10) 



where g(s), g{m 2 ) are the pion-nucleon coupling constants for pion off and on mass shell. 
This assumption can be justified by the fact that f n is weakly depending on q 2 and it is 
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in agreement with the ChPT expansion at small q 2 Q|. Therefore measuring the g n NN{ s ) 
coupling constant gives information on the axial and induced pseudoscalar FFs of the nucleon 
in the chiral limit (neglecting the pion mass). 

The matrix element is expressed in terms of the hadronic H and leptonic J currents: 

m* = — H ;j»( q ), h; = v( Pl )o; u ( P2 ), j»{ q ) = v( P+ )^u(p„), (11) 

where the index i = 0, — refers to tt° and tt~ respectively. The cross section for the case 
of unpolarized particles has a standard form (we imply that the nucleon target (proton or 
neutron) is at rest in the Laboratory frame): 

da* = — - — V \M' l \ 2 dT, p 2 = E 2 - M 2 , (12) 
16PM 1 ' v ; 

where E (P) is the energy (the modulus of the momentum), and dT is the phase space 
volume: 

1 d 3 p + d 3 p_d 3 q 7T 4 
dT = J2^^^72E- 6 iPl + P 2 ~P+-P-- ( 13 ) 

The phase space volume can be written as: 

^LdT ^ q 
2E W q (2tt) 



with 



d 3 p + d 3 p_ 4 
d^ q = n — T — 5 P-)- 



2e+ 2e_ 

Considering an experimental set-up where the pion four-momentum is fully measured, we 
can perform the integration on the phase space volume of the lepton pair: 

J dT q J2 JMMz) = ~ y(<? 2 + 2/i 2 )/5©(g 2 - V) (gr - , (14) 

where 6 is the usual step function, \x is the lepton mass and = y/l - (V/g 2 ) . 
The cross section can be expressed in the form: 

d(X = 6^ (q 2 ) 2 V 2^> (15) 

with 



s = (q n + q) 2 = 2M(M + E), r = y/l- (4M 2 /s) (16) 



and 



& = ( 9,u ~ ^f) jTr{pi - M)0;(p 2 + M)(0*)*, i = 0, -. (17) 



Using Feynman rules we can write (see Figs. dl21): 

°, = ^ (pi — q) 2 — M 2 ^ m ~ 

+ ^^9(s)F Aq 2 )^ (18) 

Note that the hadronic current J® = v(pi)0®u(p2) is conserved J®q^ = 0, but J~ = 
v(pi)0~u(p 2 ) is not conserved: 

q»J- = [(-F!{q 2 ) + F?(q 2 ))g(ml) + g(s)F w (q 2 )} v{ Pl ) l5 u(p 2 ) = Cvip^ufa). (20) 

Therefore, to provide gauge invariance, it is necessary to add to 0~ a contact term with 
the appropriate structure ([20)1 . The explicit expressions for V° and V~ are given in the 
Appendix. 

Selecting the coefficients which depend on pion energy in T>j, Eq. f!17|) . an analytical 
integration on the pion energy can be performed. In the limit of small lepton pair invariant 
mass, q 2 , after integration on pion energy, the differential cross section with respect to q 2 
becomes: 



Eq. (|2*T|) contains one of the most important results of this work, as it shows that the 
measurement of the cross section at small q 2 allows to determine experimentally the off 
mass shell pion nucleon coupling constant. 

Writing the differential cross section in the form: 

da {q 2 + 2fi 2 )f3 



<f 2 + R W) 



(22) 



dq 2 (q 2 ) 2 

with c and R(0) finite functions of s. After integration of lepton invariant mass, we find: 

a " = L P q2 = w + m s) ( iog ? - 9 + [ f ^ s) - fl(o - * (23) 

The first term in the right hand side of Eq. (|23|) is divergent for massless leptons, and 
induces a rise of the cross section (especially in the case of electron positron pair). However 
it is very hard to achieve experimentally such kinematics, q 2 — > 4fi 2 . The total cross section 
can be integrated within the experimental limits of detection of the particles. Such (partial) 
total cross section will be calculated below. 
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III. KINEMATICS 



In the Laboratory system, useful relations can be derived between the kinematical vari- 
ables, which characterize the reaction. The allowed kinematical region, at a fixed incident 
total energy s can be illustrated as a function of different useful variables. 

One can find the following relation between q 2 , the invariant mass of the lepton pair and 
the pion energy: 

q 2 = (pi +P2- q,) 2 = 2M 2 + ml + 2M{E - E n ) - 2 Pl q n = s + ml - 2E^M - 2 Pl q„, (24) 
with 

2 Pig7T = 2E W E - 2y/E%-rr%Pcoae v , (25) 

where 9 n = "p\q\ is the angle between the antiproton and the pion momenta (in the Labora- 
tory frame). 

The limit — 1 < cos^ < 1 translates into a maximal and a minimal value for the pion 
energy. The allowed kinematical region is shown in Fig. |3] left, for three values of the 
beam energy: E = 2, 7, 15 GeV. To this constraint, one should add the minimal thresholds 
q 2 > Am 2 and E n > m n . For the minimal value of q 2 ~ ml , one can plot the dependence of 
the pion energy on 8 n (Fig. El right), for different values of the beam energy. As the energy 
increases the kinematically allowed region becomes wider. At backward angles the maximum 
pion energy becomes larger at small s values. For larger values of q 2 , E n is smaller. 

For fixed values of the lepton pair invariant mass, the pion energy can take values in the 
region: 

^— = S ~ q < ±* < S ~ q = ^— (26) 
M s(l + r) ~ M ~ s(l-r) M ' V ; 

neglecting the pion mass. 

The phase space volume of the produced pion can be written (neglecting terms ~ ml/m 2 N ) 
in three (equivalent) forms: 

d^qi j 2rr 2 



dq 2 5[q 2 - 2E W (E + M - P cos 6„)]E n dE w d cos 6 n 



2irE n 

= E n dE n d cos 9k (27a) 
da 2 

= M^-dE^ (276) 

sr 

q 2 M 2 dq 2 d cos 9 w 

= hrr^ — Ai2 27c ) 

s z {l — r cos9n) 
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E n [GeV] © 

FIG. 3: (Color online) Left: The kinematical limit for q 2 is shown for cos 6 n = — 1 and for cos 6 n = 1 
in the Lab system as function of the pion energy for different values of the beam energy: £7=2 
GeV 2 (dotted line), £7=7 GeV 2 (dashed line), £7=15 GeV 2 (solid line). The allowed kinematical 
region lies below the curves. 

Right: Kinematical limit for the pion energy E n as a function of the pion angle (Lab system), for 
E=2 GeV 2 (dotted line), -£7=7 GeV 2 (dashed line), £7=15 GeV 2 (solid line), for the minimum value 
of q 2 ~ m 2 . 



IV. AXIAL AND ELECTROMAGNETIC FORM FACTORS 



Experimental measurements on FFs are object of ongoing programs in several world 
facilities. Hadron FFs are measured in space-like (SL) region through elastic electron hadron 
scattering, and in time-like (TL) region through annihilation reactions. It has been only 
recently possible to use polarization techniques. The availability of high intensity, high 
energy polarized beams allows to extend these measurements to large q 2 regions and to 
achieve very high precision. 

The theoretical effort for a complete description of the nucleon structure should be ex- 
tended to a unified picture which applies to the full kinematical region (SL and TL) jl8^ . 
Few phenomenological models, developed for the SL region can be successfully extended to 
the TL region 19]. A tentative extrapolation of a TL model to SL region has also been 
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121. 



done, and constraints have been found from the few available data 

In TL region, data for EM FFs exist over the NN threshold, up to 18 GeV 2 |20| but a 
precise separation of the electric and magnetic contributions has not yet been possible, due to 
the low statistics. Moreover, FFs are complex quantities, and polarization experiments are 
necessary to determine their relative phase. Presently only the modulus of FFs is available, 
under the hypothesis \G^\ 2 = \G^\ 2 or G E = 0. These data have been obtained in the 
reactions e + + e~ p + p (see Ref. 21 1 and references therein) and, more recently by the 



radiative return method 



22|, in the region over the kinematical threshold s > 3.52 GeV 2 . 



The possible existence of an NN resonance under threshold has also been predicted as an 
explanation of the fact that TL FFs are larger than SL ones, at corresponding \q 2 \ values. 
In order to give quantitative predictions, for the cross section of the processes (JTJ) and (j2J), 
it is necessary to know the value of FFs in the unphysical region. As data are not available, 
such estimation can be only done in the framework of models, based, for example, on VMD 
or on dispersion relations, which predict several discontinuities due to meson resonances. 
Not all nucleon models give expressions for FFs which can be extended to time-like region, 
and not all nucleon models give a satisfactory description of all four nucleon FFs. 

Following [3], in the present calculation we use two models for electromagnetic FFs: a 
model based on Ref. 2^|, which firstly predicted the behavior of the proton electric FF as 



found from recent polarization experiments [7jJ , recently extended to TL region. This model, 
as all VMD inspired models, has poles in the unphysical region, in correspondence with the 
meson resonances which are taken into account. In order to have a smooth parametrization, 
we considered also a perturbative QCD (pQCD) inspired model prescription (corrected by 



dispersion relations 



24j ) which reproduces the experimental data in TL region, does not 



show discontinuities and can be considered an 'average' expectation: 



\G N E \ = \G»\ = , -, q 2 > A 2 , (27) 

2 S_ i ~ 2 

A 



? 4 (rn^ + - 



where A = 0.3 GeV is the QCD scale parameter and A is fitted to the data. This parametriza- 
tion is taken to be the same for proton and neutron. The best fit is obtained with a parameter 
A(p)= 98 GeV 4 for the proton and A(n)= 134 GeV 4 for the neutron, which reflects the fact 
that in the TL region, neutron FFs are systematically larger than for the proton. In princi- 
ple, this parametrization holds only for very large q 2 values, but, in practice, it reproduces 
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the existing data quite well in the whole physical region. Evidently it is meaningless at small 
q 2 , (q 2 < A 2 ), and it has not the good normalization properties for q 2 —>■ 0. 

Usually the data are shown in terms of the Sachs FFs, electric Gg, and magnetic G^, 
which are related to the Pauli and Dirac FFs by the following relations: 

G N E = F» + F 2 N , G N M = F» + tF 2 n , t = q 2 /(AM 2 ). 

. They correspond in a nonrelativistic limit (or in the Breit frame) to the Fourier transform 
of the charge density (electric form factor Ge) and magnetization distribution (magnetic 
form factor Gm) of the proton. 

The behavior of these FFs is shown in Fig. compared with the existing experimental 
data. For the neutron, the first and still unique measurement in TL region was done at 
Frascati, by the collaboration FENICE Q 

Concerning the pion FF, a reasonable description exists in the kinematical region of 
interest here, for a recent discussion see Ref. [2fi| . For the sake of simplicity, we use here a 
p meson saturated monopole-like parametrization, which takes a Breit Wigner form in TL 
region: 

2 

771. 

FM 2 ) = 2 2 P ■ r • (28) 
nip — q z — ivripY p 

Data on axial FFs in TL region do not exist, and they suffer in SL region from a model 
dependent derivation. In SL region, the nucleon axial FF, GA{q 2 ), for the transition W* + 
p —>■ n (W* is the virtual P^-boson), can be described by the following simple formula 27| : 

G A (q 2 ) = G A (0)a~q 2 /rn 2 A r n (29) 

with rriA = 1-06 GeV, if n = 2. A simple analytical continuation of this prescription to 
the TL region, presents a pole in the unphysical region. Therefore we used a 'mirror' 
parametrization from SL region: 

FF^ TL \\q 2 \) = FF^ SL \\q 2 \). 

Such a prescription is, in principle, va.id only at very large f, since it obeys to asymptotic 



analytical properties of FFs |l8| . 

For comparison, another parametrization is also used 



28], inspired from [23]: 



Gf{q 2 ) = d(q 2 )G A (0) 



1 — a + a- 



rn 



(30) 
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FIG. 4: (Color online) Nucleon electromagnetic FFs in time-like region: proton electric FF (a), 
proton magnetic FF (b), neutron electric FF (c), neutron magnetic FF (d). Data are from 
and predictions from model Ql (solid line), and from pQCD (dashed line). 



The parameter a = 1.893 ± 0.02 has been fitted on the available data, in the SL region, 
averaging the dispersion due to the model dependent extraction of the data themselves, 
vtia — 1.235 GeV is the mass of a light axial meson, and d{q 2 ) = (1 — ^q 2 ) 2 is the function 
describing the internal core of the nucleon. We take 7 as a fixed parameter, from the fit of 
nucleon electromagnetic form factors : 7 ~0.25 GeV -2 . Let us note however that this value 
is not good from a t channel point of view, because it gives a pole in the physical region, 
to = ^=4 GeV™ 2 . In order to extend the expression (|30|) to TL region, following 23j, a 
phase is added to the dipole term. Moreover the complex nature of the axial FF is insured 
by adding a width to the axial meson, and we replace the propagator as in Eq. 

The expression for the TL axial FF is therefore: 



G A (q 2 ) = d(q 2 )G A (0) 



1 — a + a 



m A ~ <? 2 — i> m A^A 



(31) 
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FIG. 5: (Color online) Proton axial FF from VMD inspired model (solid line) and from a dipole 
extrapolation (dashed line). 

with d(q 2 ) = (1 -je iS q 2 y 2 , where T A =0.U0 GeV, and 5=0.925. 

The two models for Ga used in the calculation of the cross section are shown in Fig. |S] 
and differ by one order of magnitude. Moreover, en enhancement is expected from Eq. (|31|). 
in correspondance to the mass of the axial meson. 

V. RESULTS 

The differential and integrated cross sections were calculated for several values of the 
antiproton energy and different choices of FFs described above. 

The differential cross sections, Eq. (|T3jl . as a function of E n and q 2 , Eq. (27a), are shown 
in Figs. |B1 and U\ for the reactions p + p — > tt° + £ + + £~ and p + n — ► 7r~ + £ + + i~ and at 
E= 7 GeV 2 . As one can see from the figures, the differential cross sections are large and 
measurable in a wide range of the considered variables. It is reasonable to assume that the 
region up to q 2 = 7 GeV 2 , at least, will be accessible by the experiments at FAIR. 

The discontinuities in the small q 2 regions are smoothed out by the steps chosen to his- 
togram the variables. However, depending on the resolution and the reconstruction efficiency, 
it will be experimentally possible to identify the meson and nucleon resonances. 
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FIG. 6: (Color online) Left: Double differential cross section for the reaction p + p — > ir° + 1 + + £ 

fl 

as a function of q and E n , using FFs from 23] for nucleon and 1)30(1 for axial FF. Right: Same 
quantity as in the left plot for the reaction p + n — > 7r~ + ^ + + i~ . The kinematical constraints in 
the (-E/Tr, g 2 )-plane shown in Fig. EJare visible here. 




FIG. 7: (Color online) Left: Double differential cross section for the reaction p + p tt° + i + + 1~ 
as a function of q 2 and -E^, using pQCD inspired nucleon FFs and dipole axial FFs. Right: Same 
quantity as in the left plot for the reaction p + n — > 7r _ + £ + + 
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FIG. 8: (Color online) Left: Differential cross section for the process p + p — > 7r° + £ + + l~ as a 
function of q 2 , with FFs from 2j| for nucleon and (JSOJ) for axial FF (solid line) and with FFs from 
pQCD inspired nucleon FFs and dipole axial FFs (dashed line). Right: Same quantity as left, for 
the reaction p + n — > ir~ + i + + £~ . 



The differential cross section as a function of q 2 can be obtained after integrating on the 
pion energy, Eq. (fT5)l. with the help of Eq. (27b): 

da^_ a 2 f3(q 2 + 2fi 2 )M ^ 
dq 2 



(32) 



6s7rr (q 2 ) 2 sr j Ej n in 
where the integration on the hadronic term is detailed in the Appendix. The result of the 
calculation is shown in Fig. |HJ For charged pion production, the presence of the axial FF is 
the reason of a larger cross section as compared to the neutral pion case. For both reactions, 
again, the present calculation gives an integrated cross section of the order of several \xh in 
the unphysical region, for both choices of FFs. 

The q 2 dependence is driven by the choice of FFs. In case of pQCD-like FFs, the behavior 
is smooth and similar for proton and neutron. In case of FFs from j^J , the resonant behavior 
due to p, to and <fi poles appears in the figures. The partial total cross section, integrated 
over q 2 , Eq. is shown in Fig. El as a function of q m i n 3> 4p 2 for two different values 

of s, s=2 and 5 GeV 2 . Evidently the calculations for the VMD FFs 23J] have been done 
beyond the resonance region, due to the divergence of the integrals at the meson poles. 
Let us estimate the uncertainties inherent to our model assumptions. The off mass-shell 



effects were previously discussed for these particular reactions in 



has been devoted to this problem in the past (Ref. 29( and ref. therein), where it was 



] . Large theoretical effort 
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FIG. 9: (Color online) Left: Total cross section for the process p + p — > tt° + £ + + l~ as a 
function of q m in for different values of s and different FFs: dotted (solid) line for E =7(12) GeV 2 , 



nucleon FFs from 23] and axial FF from Q3U[) : dash-dotted (dashed) line for E =7(12) GeV 2 , 
pQCD inspired nucleon FFs and dipole axial FFs. Right: Same quantity as in the left plot, for the 
process p + n — ► ir~ + £ + + i~ . 



shown that indeed off-mass shell effects can be large and have the effect to increase FFs. In 
the region of virtuality p 2 < 0.5 GeV 2 it was calculated that off mass shell effects modify 
FFs at the level of 3 %. The Q 2 -dependence of FFs is not changed significantly, when one 
of the particles goes off shell. It is interesting to note that the ratio of electric and magnetic 
nucleon FF is rather insensitive to off shell effects. Moreover ChPT arguments 30] support 
a smooth behavior of FFs as a function of the degree of nucleon virtuality 5 ~ | p | < 2, 
which does not exceed 10%. For the present kinematics the virtuality involved varies in the 
interval, 2m n /M < 5 < 2 — 2/71^ jM . In the case of detection of soft pions in laboratory 
frame, errors arising from off shell effects will decrease to 3 -j- 5%. 

In the intermediate state, in principle, resonances can be excited. The estimation of 
contribution of the A resonance, compared to nucleon intermediate state, is of the order of 
MA.fi'ATVTr/ 'qnn-k — 0.07. Higher resonances have smaller contribution. 

These considerations support an estimation on the precision of our model at the level of 
10%. 



18 



VI. CONCLUSION AND PERSPECTIVES 



The differential cross section for the processes p+n — > n~ +£~ +£ + andp+p — > n +£~ +£ + 
has been calculated in the kinematical range which will be accessible in next future at 
FAIR. The main interest of these reactions is related to the possibility of measuring nucleon 
electromagnetic and axial FFs in the time-like and in the unphysical regions. As previously 
pointed out , varying the momentum of the emitted pion allows to scan the q 2 region 
of interest, keeping the beam energy fixed. 

In Ref. 2! it was also noticed that in the lepton invariant mass squared distribution, a 
divergent term was present: 

do ^ [2F^g 2 ) + F^g 2 )] 2 

dq 2 ~ (q 2 ) 2 ' q U ' {66) 

(in our notation the right hand side of Eq. corresponds to C 2 / {q 2 ) 2 )- It was argued that 
the singularity at q 2 — > cancels due to 2F^(0) = 1 and 2F 7r (0) = — 1. This compensation 
takes place if the following equality g^NN^) = 9-KNN{ m ^) holds, which is verified for anni- 
hilation at rest. In the present work, one can not rely on such assumption and the validity 
of this relation has to be verified experimentally. 

The detailed measurement of the double differential cross section, as a function of q 2 
and E n allows in principle to extract all nucleon FFs which are involved in the considered 
reactions. A precise simulation of different processes involving the production of a pion will 
be necessary with a study of the best kinematical conditions in order to minimize background 
contribution. In particular the reaction p + p —>■ 7r° + tt° has been identified as a potential 
source of background in the e + e~-spectrum due to its Dalitz decay n° — > e + e~ r y. 

The assumption about the validity of a generalized form of the Golberger-Treiman rela- 
tion, which allows to relate the pseudoscalar g^NN coupling constant to the axial nucleon 
FF can be experimentally verified in case of small invariant mass of the lepton pair. The 
possibility of measurement of heavy negatively charged pions ( tc' -radial excitation of pion, 
M' = 1300 MeV) in the process of 7r~ production can be taken into account in the following 
way. The matrix element for the excitation of a resonant state for the virtual charged pion 
in the intermediate state, M res , can be written as follows: 

M res = j s _ M g r; tMT ^ ^)^fe)% 2 ) - h (34) 
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The corresponding cross section is: 

res ~ 12m (s - M«)2 + M«r« 27rE n ' [ ' 

where A(g 2 ) is the transition from factors for the vertex -k'-k^* . We do not consider processes 
involving vector mesons, A resonances and higher excited nucleon states, estimating that 
their contribution does not exceed 10%. 

In case of multi-pion production, the quantity s\ = (pi + P2 — q) 2 — m 2 becomes positive. 
Varying si at fixed beam energy, by changing q 2 and 9 n , it is in principle possible to identify 
and study other mechanisms, as the excitation of heavy pion resonances, 7r', or the possible 
presence of a NN 'quasi-deuteron' state under the kinematical threshold for pp annihilation 
in two leptons. The study of multipion production will be the subject of a forthcoming 
publication. 
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VIII. APPENDIX 

In this appendix we give the explicit expression of the coefficients entering in the calcu- 
lation of the cross section, as well as useful integrals. 

Let us define g 2 -dependent terms, which contain FFs and the necessary constants: 

f a (s) = F v (c?)G nN s(s), f lN (q 2 ) = g(ml)F t N (q 2 ), i — 1,2, N — n,p, 
C(s)=f a (s)-f lp (q 2 )+f ln (q 2 ), 

and the quantity: 

X = MZL = ( s _ g 2 )/(2M^) - 1. (36) 
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Let us write the expressions for the hadronic part of the matrix element (|17|) : 
- for the process p + p — > £ + + £~ + n° 

^ N (1-X) 2 1 u , |2 (X+1) 2 



P° = I/: 



2pl 



E — M 1 

M 2 



g 

4M 2 



for the process n + p — > £ + + £ + 7r 



. i j,k;j<k 



X 



2\C\ 2 s 
q 2 



+ \fi P ~f: 



2p\ 



X 



(37) 



The explicit form of the coefficients is: 



Ci P ,i P - 4X, C 2pt 2 P - ( 1 + Ys X ^ ' ^ lp ' 2p 



i, j, k = Ip, 2p, In, 2n, a. (38) 



-3(1 + X), C a , a = ^-4, 
s 



1 + 



2sX 



j Cln,2n 



-3 1 + 



Ca,lp — 2, C aj 2p — ^1 — (1 + X) , Ci Pi i„ — 4, Ci P) 2n — ^— — 2X — lj 



a 



2p,2n 

n 

^aAn 



2 + 



X 2M 2 

2, C a 2n = ~ 



+ XI , C 2p i n — ^X — — — 1 

1-^4)- 



(39) 



The structure of Pj allows to select the terms which depend on the pion energy. It is 
straightforward to perform an analytical integration on the pion energy, using the following 
integrals: 

(40) 



I/lli <!.r 

t 

f) It'll M 



™ ax dE n r(s - q 2 ) . , s - g 2 
— ro; o - 



2M 2 



2M 2 



X 



£ - dE n 1 _s- q 2 
E min ~WX ~ 2.M- 



. 1 + r 

in r 

1 — r 



6(*-r); 



(41) 



with r given in (jlfij) and £ = ln[(l + r)/(l — r)]. The result of the integration on the pion 
energy is: 

- for the process p + p — > £ + + £~ + n° 



rpmax , _ 

•pi) u / - ■ 



M 



- = 6 2|/ lp -/ 2p | 2 £+|/- 



2p| 



E-M 



M 



-r + 1 



4M 2 



(2r - £) 



for the process p + p — > £ + + £ +7r° 



P~ 



M 4 



E 2 • 2 E KitMfiK) + \c\ 



,2rs 



(42) 



(43) 
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where 



K 

ly In, In 



K, 



a,lp 



K- 



2p,2n 



-3£, K a , a = [ ^- - 4 ] r, 
4(£-r), tf 2 „ 2n = ^(r + £(£-r)), K ln , 2n 



-3£ 



2r, K ai2p 



9 

!- — )£, = 4r, i^i p ,2„ 



3£- 
-2r, i^ a o n 



2M 2 / 



K- 



2p,ln 



1-21 



(44) 
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